On the nucleon self-energy in nuclear matter 
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. We consider the nucleon self-energy in nuclear matter in the absence of Pauli blocking. It is 

' evaluated using the partial-wave analysis of AW scattering data. Our results are compared with 

, that of a realistic calculation to estimate the effect of this blocking. It is also possible to use our 

results as a check on the realistic calculations. 
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A fundamental problem of nuclear physics is to explain the properties of nuclear matter (and finite nuclei) in terms 
of an effective field theory at low energy based only on the (chiral) symmetry of QCD. While such a theory has been 
qJq ' eminently successful for systems like tttt and ttN [1], a satisfactory theory for the NN system has been difficult to 
^f) , formulate due to the presence of two- nucleon bound or virtual states close to the threshold of NN scattering [2]. 
' In particular, the leading chiral four-nucleon interaction predicts an absurdly large value for the self-energy of the 
, nucleon at normal nuclear density [3]. 

i There is, however, a semi-phenomenological approach that yields fairly accurate values for different observables 
in nuclear matter. Here the NN potential is constructed by exchanging low mass bosons in the t-channel [4]. The 
coupling and other parameters in the potential are determined by experimental data on the deuteron and the low 
G ■ energy AW scattering. The dynamics is formulated on the basis of a relativistic version of the Brueckner-Hartree-Fock 
I , method [5,6], where the reaction matrix satisfies a three-dimensionally reduced Bethe-Salpeter equation in nuclear 
^ ' medium. The Dirac equation for the nucleon incorporates the scalar part of the self-energy due to its interaction with 
^ nucleons in the medium. The self-energy itself is given by the diagonal element of the reaction matrix. The system 
^ of equations is then solved self-consistently. 

^ ' In this work we study the nucleon self-energy in a certain theoretical limit. We observe that if we suspend the Pauli 
blocking operator (projecting onto the unoccupied states) in the equation for the reaction matrix, it coincides with 
the one for the scattering matrix in vacuum. Further, if we do not include the relevant part of the single particle 
self-energy in the mass term in the Dirac equation, the requirement of self-consistency does not arise any more. In 
this limit the self-energy is given by an integral over the spin-averaged, forward AW scattering amplitude in vacuum, 
which can be evaluated entirely with the experimental data. 

There is a well-known expansion in statistical mechanics, called the virial expansion, whose first term for the in- 
medium self-energy would give exactly the theoretical limit considered above [7-9]. Here we employ this method to 
derive the formula for the limiting self-energy of the nucleon. We then evaluate it at different nuclear densities, using 
the phase-shift analysis of WW scattering data, independently of any WW potential. It is then compared with the 
realistic calculation [6] to assess the importance of the effect of Pauli blocking in nuclear matter. 

Our calculation would also serve as an important check on the realistic calculation. One has just to repeat the 
calculation of the self-energy in the original framework itself, using the phenomenological potential and the physical 
nucleon mass, but in the theoretical limit of omitting the Pauli blocking operator in the reaction matrix. This result 
as a function of the nuclear density must agree with that of the present calculation. This, in turn, would confirm our 
assertion that it is indeed the Pauli blocking effect which distinguishes the realistic calculation from the one presented 
here. 
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II. DERIVATION OF SELF-ENERGY FORMULA 



Here we obtain the leading term in the virial expansion for the nucleon self-energy in nuclear medium [7-9]. We 

begin by stating clearly the normalization of different quantities. Omitting the nucleon isospin index [10], wc take the 
creation and the annihilation operators for the nucleon with momentum p and spin projection u (= ±i) to satisfy 
the anticommutation relation, 

{bip,a),b^ip',<j')} = {2TTf2EpS{p-p')6„„,, Ep = ^p^+m'. (2.1) 

The single nucleon state is defined as |p, (t) = cr)|0) . Similarly the two nucleon state is [p^, cti; ^2) = 

h'^{Pi,ai)b\p2-,(T2)\^)- Clearly their normalization is fixed by the anticommutation rule (2.1). The (positive energy) 
Dirac spinors are normalized such that the spin sum over these spinors, to be used below, is given by 

^ u{p, a)u{p, a)='i) + m . 

(7 

The derivation starts by considering the nucleon self-energy in vacuum, which may be expressed as an S'-matrix 
element, 

-i{2w)^S\p[-p,)u{p[,a[)^^°\pMp„a,) = mp[,a[) {S - l)b'<{p„a,)\0) , (2.2) 

where S is the familiar scattering matrix operator, S = Te^-f '^''^^^^^'^^^ , for an interaction Lagrangian Cmt- The 
subscript 1 on the variables of the particle anticipates another one with which it will interact in the medium. In fact, 
we shall express below the nucleon self-energy in nuclear medium in terms of the NN scattering amplitude in vacuum, 
defined as usual by 

{p'i,'^[;P2,'^2\S - l\Pi,(Ti;p2,cT2) = ii2TT)'^S'^{p[ + P2 - Pi - P2) M{pi,c7i; p2,(T2 p[,a[,;p'2,a'2) , (2.3) 

where M stands for the scattering matrix sandwiched between spinors corresponding to the final and the initial states 
of the two nucleons. 

Below we shall meet the spin averaged amplitude in the forward direction, 

M{pi,p2^Pi,P2)) = ^'^ M{pi,ai;p2,(T2 ^ Pi,(Ti;p2,(T2), (2.4) 

(Tl,<T2 

With our normalization of states, the amplitude M is Lorcntz invariant. 

To obtain the self-energy in nuclear medium, we have to replace the vacuum expectation value in Eq. (2.2) by an 
appropriate one. Although we specialize later to zero temperature, we take here the most general average over an 
ensemble of systems maintained at temperature T{= 1//3) with nucleon chemical potential fi. Thus the in-medium 
self-energy E is given by 

^z{27T)^6\p[-pi)u{p[,a[)npHPi,ai) = {b{p[, a[){S - l)b{p„ a^)) , (2.5) 

where for any operator O, 

Here H is the Hamiltonian and Af the nucleon number density operator. Clearly this form of the Boltzmann weight 
breaks explicit Lorentz invariance and singles out the rest frame of the medium [11]. 

We now make use of the virial expansion to first order in density. For an operator O, the ensemble average in 
nuclear medium can be expanded as 

(O) = {0\O\0)+J2 J j^^n{E,,)(j>2,<J2\0\p2,cJ2) , 

where n{Ep) is the nucleon distribution function, n{Ep) = l/[e'^(^f -|- 1]. Applying it to the left hand side of 
Eq. (2.5), we get for the difference S(")(p) = S(p) - ^'^^Hp), 

-i{27^)'5\p[-pi)u{p[,a[)J:^"HPiMPi,<ri) = Yl I (^^^5^n(£^pJ(P2,'^2|Kp'i,'^'i)(5-l)&^(Pi,'^i)|P2,'^2) . 

(2.6) 
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The matrix clement in Eq. (2.6) will be immediately recognised to be the AW scattering amplitude defined above by 
Eq. (2.3). Cancelling the ^-function on both sides, we set a[ = ai and sum over ai also to get 

-tr{S(")(pi)(^ +m)} = 4j .^-f^n{E,,)Mipi,p2 - pi,P2), (2.7) 

where the tr(ace) is over matrices in Dirac space. Note the similarity of this equation with the corresponding one 
in Brueckner theory [12]. There is, however, one important difference: Our first order formula has the scattering 
amplitude in vacuum,, while it is the amplitude in medium that enters the equation in Brueckner theory. We shall 
discuss this point again in Sec. IV. 

So far we did not state explicitly the isospin structure of the amplitude M, which is now easy to figure out. We 
consider symmetric nuclear matter and work in the limit of iso-spin symmetry. Let the traversing nucleon be in any 
one of its isospin states, say a proton. It may scatter with a proton or a neutron in the medium. The amplitude is 
therefore given by the sum, 

M = Mpp^pp + Mpn^pn ■ (2.8) 

We now restrict to the case, where the three-momentum is set equal to zero. Then the rest frame of the medium 
coincides with the lab frame of the scattering process. The self-energy in this frame has the simple Dirac matrix 
structure, 

S(") = U-1 + V'y°, (2.9) 

where the coefficients U and V depend only on the nucleon density. Then the left hand side of Eq. (2.7) simplifies to 
— 4m(?7 + V). On the other hand, the nucleon propagator with self-energy correction, i/{]l>i — m — Y,^"'\pi)}, reduces, 
for Pi — 0, to 

\l+7°), (2.10) 



Po-{m + U + V)2 
in the vicinity of the pole. The shifted pole position is thus given by 

m*-'j = m + U + V = m-^ f '^'f' n{p2)M(p2), (2.11) 

where m* is the effective mass of the nucleon and 7 gives the damping rate of nucleonic excitations. 

III. EVALUATION 

In our evaluation we restrict ourselves to nuclear matter at zero temperature. In this limit we assume the nuclear 
medium to be a non- interacting Fermi gas with all states filled upto the Fermi momentum pp, so that n{p) 
^{Pf — \p\)- For the symmetric medium the number density is then given by 

-MpF - \p\) = 



where is related to the chemical potential by = \/ \^ — rm?. 

The scattering amplitudes are generally analysed in the centcr-of-mass (cm.) frame, where they arc normalized in 
a slightly different way to get a simple expression for the differential cross section. One defines a scattering amplitude 
/ related to M by 



da _ |M|2 



dQ. {SttE) 



= i/r (3.1) 



where E is the total energy in the cm. frame. The results of partial wave anal ysis are gener ally given as functions 

of the lab kinetic energy T{= \Jp\^ m? — m), in terms of which we have E = ^J2m{2m + T). We also note here the 
expression for the cm. momentum in terms oi T, k = y' mT/2. 
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We can now write Eqs. (2.11) as 



m 



7 = TO--W— / dTVT(2m + T)f(T) , (3.2) 

TT V TO Jo 



where Tp, the upper Umit of the integral, is related to pp by pp = ^JTp{2m + Tp) and 

f / pp — >pp ^ f pn — >pn ' 

the bar indicating spin averaging as in Eq. (2.4). In terms of amplitudes with definite iso-spin, 

7 = 3/2/^=^^l/2/^=°\ 

The full scattering amplitudes are expanded in a series of partial wave amplitudes, which may then be determined 
by fitting with experimental scattering data. For the spin averaged, forward scattering amplitudes, this expansion 
takes a particularly simple form [13,14], 

-f\E) = 2.\Y.i2j + l)fli\E) (3.3) 

jsl 

Here the factor of 2 takes into account the identity of the two nucleons in the scattering process. The total angular 
momentum j is obtained by coupling the total spin and orbital angular momenta s and I respectively. The Pauli 
principle restricts the possible amplitudes by requiring the quantum numbers to satisfy 

(_i);(_i)i-(_i)i-/ = _i 

In general, the amplitude /^■'"(-E) is a matrix in the I space, whose diagonal elements enter the sum in Eq. (3.3). 
Below we shall remove the superscripts / and s on the partial wave amplitudes. 

The form of the partial wave amplitudes is determined by the unitarity of the S'-matrix. Thus for the uncoupled 
waves, s = Q,l = j and s = 1, / = j, we have the single element. 



fl^{E) ^ fl{E) = (e^^^- (i^) _ i) |2^k, l = j, 



where 5j is the phase shift, a real function of E. But the waves s = 1,1 = j ±1 are coupled, leading to a 2 x 2 matrix 
amplitude with the diagonal elements of the form [15], 



fl^iE) = (e2^^?(^)cose,(£) - l) /2ik, l=j±l. 



where we have the mixing parameters ej {E) in addition to the phase shifts. 

We may now evaluate the integral (3.2), using the phase shift analysis of the Nijmegen group [16]. Alternatively 

we may take advantage of the reconstruction of the full (Saclay) amplitudes from this analysis, also carried out by 
the same group. The general amplitude may be written in the cm. frame as a 4 x 4 matrix in the Pauli basis as [17] 

■M.{p,p') = ^{as + bs + {tts - bs)(Ti-ncr2-n + {cg + ds)cri m<T2 m + {cg - ds)a-i-lcr2-l} (3.4) 

where the Saclay amplitudes, a^, bg , Cg and dg, are complex functions of the energy and scattering angle. (We omit a 
fifth amplitude, which is zero in the forward direction.) Here I , m and n arc three mutually orthogonal unit vectors. 
The Pauli matrices <ri , <T2 act on the Pauli spinor x's of the first and second nucleon. The spin-averaged, forward 
amplitude / is obtained from M. as 

1{E) = i ^ xUxlMMXa.X.. = \{as{E) + bg{E)) (3.5) 

Cl,(T2 

With the values of the Saclay amplitudes [16], we evaluate the real and the imaginary parts of the integral (3.2) at 
different Fermi momenta. The results are shown by the solid curves in Figs. 1 and 2. 
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FIG. 1. Shift in nucleon mass in nuclear matter as a function of Fermi momentum. The sohd curve results from the 
partial-wave analysis of the Nijmegen group, while the dashed one is from the s-waves in the effective range approximation. 




FIG. 2. Damping rate of nucleonic excitation in nuclear matter as a function of Fermi momentum. The origin of solid and 
dashed curves are the same as in Fig.l. 

For an independent, but approximate, estimate, we also evaluate the integrals by including only the s-waves in the 
effective range approximation. Here an s-wave amplitude is written as /o = cotS — ik) with k cotS = ~a~^ +rk'^ /2, 
where a and r are the scattering length and the effective range. The values of these constants are long known [13]; 
For the spin singlet state, a = —23.7, r = 2.7 and for the spin triplet state, a = 5.39, r = 1.70, all in units of fm. 
This evaluation is shown by the dotted curves in Figs.l and 2. It is seen that the higher partial waves contribute little 
up to about pf = l/m~^. 



IV. DISCUSSION 



Here we have considered the nucleon self-energy in nuclear matter, in the limit of ignoring the effect of Pauli blocking 
on it. This self-energy can be expressed in terms of the forward spin-averaged AW scattering amplitude in vacuum. 
We calciilatc its real and imaginary parts, using the phase- shift analysis of experimental data on AW scattering. 

Our results may be compared with that of the self-consistent Hartree-Fock calclation [6] to get an idea of the 
importance of Pauli blocking in the Fermi medium. (Since we do not include the relevant part of the self-energy in the 
nucleon mass, it would be appropriate to compare our results with their so-called 'non-relativistic' version of results.) 



In both the calculations the mass-shift is a strongly dependent function of the nuclear density. At normal nuclear 
density {pp = 1.35 fmr^) they find the (real part of the) mass shift to be —87 MeV; in our calculation this value is 
attained at a higher density corresponding to pp = 1-70 fm~^ . 

Our calculation, which is based only on experimental data, readily provides a check on the the original relativistic 
Brueckner calculation [6]. We just need redo this calculation without the Pauli operator in the equation for the 
reaction matrix. The resulting calculation with the phenomenological potential should yield the same functional 
dependence of the self-energy on nuclear density as we find here. 
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